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There is growing observational evidence for a population of supermassive black holes (SMBHs) 
in galactic bulges. We examine in detail the conditions under which these black holes must have 
originated from primordial black holes (PBHs) . We consider the merging and accretion history expe- 
rienced by SMBHs to find that, whereas it is possible that they were formed by purely astrophysical 
processes, this is unlikely and most probably a populations of primordial progenitors is necessary. 
We identify the mass distribution and comoving density of this population and then propose a cos- 
mological scenario producing PBHs with the right properties. Although this is not essential we 
consider PBHs produced at the end of a period of inflation with a blue spectrum of fluctuations. 
We constrain the value of the spectral tilt in order to obtain the required PBH comoving density. 
We then assume that PBHs grow by accreting quintessence showing that their mass scales like the 
horizon mass while the quintessence field itself is scaling. We find that if scaling is broken just before 
nucleosynthesis (as is the case with some attractive non-minimally coupled models) we obtain the 
appropriate PBH mass distribution. Hawking evaporation is negligible in most cases, but we also 
discuss situations in which the interplay of accretion and evaporation is relevant. 



I. INTRODUCTION 

There is growing evidence for the presence of super- 
massive black holes (SMBHs) in the center of most galax- 
ies (l]-§ including our own B (but see H for a more 
skeptical view) . The origin of these black holes (and their 
relation to the host galaxy) is far from certain but sev- 
eral theories have been advanced. Of relevance is the 
the observational fact that there is a proportionality re- 
lation between the mass of the nuclear black hole, Mf,/,,, 
and that of the bulge |^J^. The bulge mass appears to 
be about a thousand times larger the black hole mass, a 
relation which holds over 3 to 4 orders of magnitude. 

Nevertheless this close relationship may be interpreted 
variously, and ultimately one is confronted with a chicken 
and egg problem - what came first: the host galaxies or 
the SMBHs? It is not inconceivable that SMBHs are 
purely the result of the internal galactic dynamics and 
their merging history; yet no one has proposed a con- 
crete mechanism for converting stellar mass objects into 
objects 6 to 10 orders of magnitude larger (with the pos- 
sible exception of ||l^). But it could also be that cen- 
tral black holes preceded any luminous activity and that 
black holes led the formation of the first galactic bulges 
and Quasar (QSO) activity |^,| If the latter is true 
one must then find an explanation for the origin of the 
primordial population of black holes. 

In both scenarios it is inescapable that black hole 
growth has taken place in the recent cosmic history. Even 
if the budge matter is well virialized, and whether or not 
it fuels the SMBH, every time galaxies merge and their 
nuclear black holes coalesce part of the bulge matter ends 
up in the central black holes Thus, as galactic merg- 
ing proceeds, the comoving density of SMBHs decreases 
and their masses increase. In Section || we spell out the 
uncertainties of this chicken and egg process, identifying 
the conditions under which a primordial population of 



black holes is necessary. We then compute the density 
and mass profile required of this population, in order to 
explain the observed SMBHs. 

The rest of our paper is devoted to proposing a cosmo- 
logical mechanics for producing the required pre-galactic 
black holes. According to our theory SMBHs are descen- 
dants of primordial black holes (PBHs) produced in the 
very early universe. PBHs are produced, for example, at 
the end of inflation |16|, in double inflation scenarios Jl7[ , 
or in first order phase transitions |18|. To fix ideas, and 
although this is not essential for our paper, we shall fol- 
low |]16[ and assume that PBHs are produced at the end 
of a period of inflation with a blue spectrum of fluctua- 
tions (with the possibility of a running spectral index Us ). 
Whatever their origin, all PBHs previously considered in 
the literature are much lighter than SMBHs, with masses 
of the order M « M0(T/ IGev)-^ for PBHs formed at 
temperature T. Hence only a very unrealistic phase tran- 
sition, at T « IMev could produce SMBHs with masses 
of the right order of magnitude. 

However the standard argument assumes that for all 
relevant cases PBH masses remain constant once they're 
formed, since evaporation and accretion can be ne- 
glected p9|-^. We show that this is not necessar- 
ily true and focus on scenarios in which significant ac- 
cretion occurs during the lifetime of PBHs. Specif- 
ically we assume that the universe is pervaded by a 
quintessence field p^-p8[|. Black holes cannot support 
static scalar fields in their vicinity and will try to "eat" 
them; quintessence is no exception. In the process their 
mass increases, so that the seeds which led to the SMBHs 
we observe today could be PBHs which have eaten too 
much quintessence. 

We examine this possibility in Sections HI to ^ In 
Section HI we estimate the effects of evaporation and 



accretion in the presence of a quintessence field. In Sec- 



tion [V we compute the comoving density of PBHs in our 
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model. Finally in Sec. |^ we compute the PBH mass spec- 
trum, and constrain the free parameters in our model in 
order to fit the requirement derived in Section ||. 



II. GALACTIC BLACK HOLE ACCRETION AND 
MERGER HISTORIES 

As outlined above, the close relationship between the 
mass of the central galactic black hole M^h and that of 
the galactic bulge may be interpreted variously. It could 
signal that the central black hole was formed by the in- 
flow of bulge matter (stars and gas), but it could also be 
that the central black hole was there initially and 
led the formation of galactic bulges. Most probably there 
was a combination of the two scenarios, and the central 
black hole could indeed be primordial, but at the same 
time it was also fed by outside matter. 

In this section we show that such a combination is in- 
deed highly feasible. We devolve the merger history of 
the galactic halos [^,0, and use a simple merging/ ac- 
cretion prescription for the behaviour of the central black 
holes iQ. We thereby show that the majority of galac- 
tic black holes today could have originated prior to the 
formation of significant halos, whilst still being in agree- 
ment with the observed mass evolution seen in QSO's 

To devolve the galactic black holes, we use a merger 
tree approach ||2^,^ to establish the histories of the ha- 
los and the supermassive black holes within them. The 
method involves prescribing a number density distribu- 
tion of halos nowadays using the Press-Schecter formula 
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An adaptation of the basic equation provides the condi- 
tional probability that a halo of mass Mg at redshift zq 
evolved from a progenitor of mass Mi = Mq — AM at 
redshift zi — zq + Az 29 1 , 
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We use the notation of Lacey and Cole where in equations 
(H) and (H), p{z) is the energy density, a'^{M) is the vari- 
ance of density fluctuations on a spherical scale enclosing 
a mass M, and u}{z) = 6c{l + z) where 5c ^ 1.686 is the 
over density threshold for density fluctuations to collapse. 
We take uj{z) from a modelled matter power spectrum us- 
ing Hq = Ibkms'^Mpc'^ and flc = 0.25, fib = 0.05, and 
il^ = 0.7. The redshift step Az is mass dependent and 



represents a realistic merging timescale for the chosen 
halo, we follow and take Az = Auj{M)/6c — 1 where 
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As a halo is deconstructed only those with a mass 
greater than a limiting mass Mi are traced, any pro- 
genitor with M < Ml is treated as unbound matter 
accreted onto the halo. We consider a range of limit- 
ing mass scales Mi = 10^ - 10^°Mq. This is akin to 
putting in a lower threshold on the halo's velocity dis- 
persion cr* ~ Vyir = {GM/ryiry/^ of ^ 40 — 50kms~^ 
where ryir is the virial radius and where we assume a 
spherical bulge for the halo so that M = 47rrj5j^p^a;o/3. 
We initially restrict our analysis to typical galaxy mass 
scales today M ~ 10^° - IO^^Mq- 

We assume that a central black hole is present in 
all halos M > Mi. To assign black hole masses in 
each halo today, we use the powerful correlation re- 
cently found between the black hole mass, Mf,/,,, at the 
center of galaxies and the line of sight velocity disper- 
sion cr* [|l^,^. Specifically we use the best fit relation, 
Mbh = 1-2 X lO*^M0(cr*/2OOA:TOS-i)3-^5 g 

The evolution of black holes through merging events 
and accretion of halo matter is a complex one, the 
timescales over which merging occurs will be intricately 
dependent on, amongst others, the size of halos and the 
ferocity of the merging event. While the mechanism for 
accretion will be dependent upon the halos properties 
(for example redshift and halo velocity dispersion). We 
do not endeavour to model these complex processes here 
and instead take a simplistic approach focused on current 
observational constraints. We assume halo mergers are 
violent events allowing black holes to merge on timescales 
significantly less than the time between halo mergers. We 
then model accretion using a relation proposed by p^ ] 
whereby a fraction of the halo gas is accreted onto the 
black hole. 

Mace = eacc(l + z)'^ Mhalo cxp h (d, /SOOA^ms"!)*] (4) 

where eacc is the accretion efficiency factor, and the ve- 
locity dispersion dependence is introduced to account for 
the reduction in the halo gas's ability to cool in the grav- 
itational well around the centre at late times. 

We are interested in the prospect that primordial black 
holes are progenitors for the SMBHs nowadays. We as- 
sume that only halos with M > Mi contribute signifi- 
cantly to accretion onto the central black hole. Subse- 
quently, if a halo only has progenitors of M < M; and 
retains a black hole then the black hole is assumed to 
have not undergone any gas accretion at higher redshifts 
and is treated as primordial. 

Figure ^ demonstrates, for the case of a single halo, 
how the accretion efficiency, eacc, in @ has a strong influ- 
ence on how much of the black hole mass could be present 
at higher redshifts. In figure we show the evolution of 
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FIG. 1. Evolution to give the present day black hole in 
a 5 X 10^^ Mq halo with redshift for various accretion efh- 
ciency factors. From top to bottom the evolutions are for 
Eacc = 10~^3 X 10-^10~^3 x 10~■^and 10"^ 




FIG. 2. Evolution of mass and number density for galac- 
tic black holes for 4 different accretion efficiencies. The data 
points are the inferred mass density values of Chokshi and 
Turner. The initial mass density and number density of black 
holes are taken to be 1 x lO^MgA/pc"^ and 4 x lO'^Mpc'^ at 
z = respectively. The number density at early time reaches 
9 x Mpc~^ at z = 10 for the scenario that is consistent 
with observations. As one views the number density evolu- 
tion from high z, the increase in number density arises from 
the formation of astrophysical black holes due to halo matter 
accretion. The contribution of astrophysical black holes be- 
comes more prevalent as the accretion efficiency is increased. 



the total comoving energy and number densities for black 
holes with redshift. There is an intuitive play off, shown 
in figure ||, between the accretion efficiency and the frac- 
tion of black holes originating at early times. The higher 
the accretion rate, the higher the chance is of black holes 
being formed at late times, during the halo merging, as 
opposed to being primordial. As discussed in p^ ], con- 
straints can be placed on the efficiency coefficient using 
inferred accretion rates from QSO luminosity functions 
p2t . One can see that the QSO evolution data places 
a tight constraint on the accretion efficiency parameter: 
for Ml = 101°Mq one requires Cacc ^ 1-85 x 10"^. In this 
scenario we see that the vast majority of black holes, us- 
ing this accretion prescription, are present prior to halo 
merging. 

An alternative scenario extends the galaxies able to 
support black holes down to dwarf scales. This can be 
modeled by lowering the mass threshold to 10^ Mq. In- 
evitably this will alter the evolutionary profile; accretion 
in smaller galaxies will be comparatively more efficient 
under this model as their gravitational well is not as pro- 
hibitive. As is shown in figure ^ this leads to a slight 
reduction in the predicted comoving energy density of 
black holes at early times. However, the general conclu- 
sions are the same in that the PBHs still contribute a sig- 
nificant fraction of the energy density today and comprise 
the majority of the total number of black holes involved 
in the evolution. 

In figure ^ we provide the main output of this section. 
We show the expected mass distribution for the primor- 
dial black hole for Mi = 10^^ Mq and Mi = lO^Afg. 
Their masses, typically Mth ~ 10"^ — lO^M©, are still 
many orders of magnitude larger than stellar mass scales, 
even considering the potential merging of early massive 
stars. We propose that such black hole masses could be 
generated through early time accretion of a quintessen- 
tial scalar field onto PBHs. In the following sections we 
investigate such a scenario and show that mass distribu- 
tions such as those is figure ^ can be obtained. 



III. THE EFFECTS OF EVAPORATION AND 
ACCRETION ON PRIMORDIAL BLACK HOLES 

There has been intense discourse regarding whether or 
not PBHs are capable of accreting radiation. Carr and 
Hawking gave arguments for negligible accretion p9[ | , but 
these were later disputed (see pO| and references therein) . 
In any case these arguments only apply to perfect fluids 
and a scalar field (j> is not a perfect fluid (even though 
an isotropic, homogeneous scalar flcld does behaves like 
a perfect fluid). In Appendix 2 we show that even for the 
simplest potentials V{4>) the field is indeed absorbed by 
the black hole leading to a mass increase with rate: 



^ = An{2M)m 



(5) 
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FIG. 3. Evolution of mass and number density of black 
holes in galaxy scale halos for 3 different limiting mass thresh- 
olds. Reading top to bottom for mass density and bottom to 
top for number density Mi = 10^°, 5 x 10^ and IO^A/q with 
accretion efficiencies of e^cc = 1.85, 1.63 and 1.38 x 10"" re- 
spectively chosen to be in agreement with low z QSO obser- 
vations 



log(M,{Mj) 

FIG. 4. Distribution of black hole masses at redshifts pre- 
ceding halo merger activity for the models in fig 3 with 
Ml = 1O^°M0 (fuU line) and 10^ (dashed line). The dis- 
tributions are normalised to their peak values. 



where (pc is the cosmological solution for (j). For more 
general potentials it may well be the case that the pro- 
portionality constant in the equation receives corrections 
of order 1, but these shall not matter for the rest of the 
argument in this paper. 

Prom eqn. ^ we note that kinetic (but not potential) 
scalar energy leads to black hole growth. This is consis- 
tent with the result that the presence of a cosmological 
constant does not lead to equivalent growth. 

Assuming now a potential of the form V = Vqc^^"^ , we 
have 
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leading to the asymptotic mass: 
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for black holes with a mass smaller than a critical mass 
Merit — t/K. These black holes eventually stop accreting 
and therefore are subdominant with respect to all others. 
For M = Merit BHs grow like t. Above this value the 
black holes would seem to grow faster that t, however 
clearly the approximations used must break down at this 
stage. In 1 20 it was shown that causality constrains these 
to grow like t as well, a result we shall use for the rest of 
our calculations. 

Hence there is a critical mass at time of formation when 
PBHs may grow proportionally to the horizon mass. This 
critical mass separates those black holes which will be of 
relevance for out scenario and those which won't. 

In addition PBHs may experience significant evapora- 
tion, via Hawking radiation. This leads to a decrease in 
their mass, at rate: 
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This equation can be integrated into 

M= {Mo - 3a{t ~ to))^^^ 
implying an evaporation time of 



(11) 
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FIG. 5. The evolution of black holes ignoring evaporation. 
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For black holes formed at temperatures much smaller 
than the Planck temperature this effect can be ignored. 

Considering now both accretion and evaporation the 
black hole mass rate equation becomes 



dM 



a 

'M2 
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(16) 



which does not have analytical solutions. However, if 
the temperature at which the black holes are produced 
is not close to the Planck temperature the interplay of 
accretion and evaporation is very simple. Black holes 
with Mo > Merit will grow with M oc i, and since their 
mass was never too small they never experience signifi- 
cant evaporation. Black holes with Mq < Merit will stop 
growing at some Moo, following (^). If this is smaller 
than lO^^g they will evaporate before today following 
( p^ ) . Accretion and evaporation happen at very distinct 
times, so although we do not have an exact solution to 
( p^ ) it is an excellent approximation simply to glue to- 
gether back to back (^) and ([l^). 

For reheating temperatures of the order of the Planck 
temperature the situation is more complicated as evapo- 
ration may be significant while black holes are accreting. 
The result of a numerical integration is plotted in Fig. g 
(to be compared with Fig. ||). We see that the effect 
of evaporation is then to shift upwards the critical mass 
Merit above which the black hole mass scales like t. In 
addition, black holes with masses lower than Merit are 



FIG. 6. The evolution of black holes formed at Planck time 
with a = a — 1. 



quickly removed by the effects of evaporation even while 
they are accreting. 



IV. PBH FORMATION AND THEIR COMOVING 
DENSITY 

Having identified the conditions required from PBHs 
in terms of early-time accretion and late-time merging, 
we now proceed to construct a cosmological scenario in 
which they are satisfied. Although not compulsory, for 
definiteness we consider an inflationary scenario with a 
reheating temperature Tr and a tilted spectrum of scalar 
perturbations with a spectral index Ug > 1. The index 
Us need not be a constant and indeed many inflation- 
ary models predict a running tilt, varying from scale to 
scale. Bearing this in mind we stress that the constraints 
upon Us discussed in this section refer to very small scales 
(the horizon size at Tr), widely different from the scales 
probed by CMB fluctuations or large scale structure sur- 
veys. 

As shown in [ p^ for suitable Us one may have produc- 
tion of PBHs in a short time window immediately after 
reheating. Typically a black hole is formed if the den- 
sity contrast on the horizon scale exceeds a given critical 
value, S > 5c- Its mass is given by: 



Mbh^kMniS-Sey/''. 



(17) 



Numerical studies with PBHs formed in a pure radiation 
background have identified 5e, k and 7 [ p2|j23| ; it may be 
that quintessence modifies these values slightly however 
for simplicity we take the calculated values k = 3.3, 7 = 
0.34. But given the uncertainty we consider two values of 
5e in this section 6e = 0.25, 1. The correct value should 
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be somewhere in between. In section ^ we assume the 
value 5c = 0.67 to outhne a precise example. 

Note that the horizon scale is important because it is 
the relevant Jeans scale for radiation but also for the 
quintessence field. In Appendix 1 we present a simple 
model for the formation of PBHs with quintessence: some 
peculiarities are found, but they should not affect the rest 
of the argument. 

The mass variance on the horizon scale at temperature 
Tis 0: 



VMev / 



(18) 



To a good approximation, as shown in figure ^ we 
may assume that all PBHs are formed immediately after 
reheating, since an then decreases sharply making black 
holes rarer. 

We have seen that we expect no black hole coalescence 
for redshifts higher than z « 10. Hence the comoving 
density of PBHs can be directly related to the probability 
of an accreting black hole being formed, Pacc- This is the 
probability oi 5 > Scrit, i-G- the probability for black hole 
formation and accretion in each horizon. 
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(where the last approximation comes from [p3[) and 
where Scrit is the value of ^, in equation (pT[), that would 
create a black hole of critical mass, Merit- The comoving 
density of accreting black holes is n « Pacc/ {^T^r^) where 
th is the comoving horizon radius at Tr- The latter is 
given by rn = 2tao/a « 80(T/Mev)-iPc so that 



71 = 4.6 X 10' 



Pa 



(20) 



Combining ([l8|), ( p^ ) and (po|), and to be consistent with 
Fig. requiring that n w 9 x 10~'^Mpc~'^ we obtain a 
value for Ug for each value of Tr. In Fig. |^ we plot the 
required ns{T) for the two values 5c = 0.25,1, using a 
fixed value of McritW- As mentioned above, the correct 
value of 5c should lie somewhere in this region. In Fig. |^ 
we plot the required tilt for three different critical masses 
with fixed 5c = 0.67. 

The tilts required fit within the constraints of , and 
can accommodate the recent reionisation tilt constraint 
of 1^ with Us < 1.27 with a high reheat temperature 
>- lO^^Ge^;). 

Note that the primordial black hole comoving density 
imposed in our considerations is much smaller than that 
considered in Indeed the black holes studied b y |1^ ] 
do not grow (as opposed to ours ). The idea in [^| is 
to produce a larger density of much lighter black holes, 
suitable to promote them to candidates for dark matter. 
Our purpose is to generate a much lower density of much 
heavier PBHs, so that they could supply the primordial 
eggs for the merging history of SMBHs. 




log(T/Gev) 

FIG. 7. The value of the spectral index Ua (on the scale of 
the horizon at the time of PBH formation) resulting in the 
correct comoving density for SMBH seeds. The two curves 
correspond to 5c=0.25,l and A = 3 is taken as an example. 
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Log(T/GeV) 

FIG. 8. The value of the spectral index Us (on the scale of 
the horizon at the time of PBH formation) for three different 
critical masses relating to A = 2, 3, 4. 5c is fixed =0.67. 



6 




log(a) [ag-l 
FIG. 9. The demography of black hole formation with a 
reheat temperature of Tr = IQ^^Gev and tilt Us — 1.33. The 
crosses represent Monte Carlo simulated black holes. The full 
and dashed lines are the horizon and critical mass, the latter 
assumes A = 2. 



V. THE PBH MASS SPECTRUM 

As an example we calculate the mass spectrum follow- 
ing the methods of to evaluate the initial mass spec- 
trum. This involves using a Monte-Carlo technique to 
obtain black hole masses which, unfortunately, is compu- 
tationally unfeasible for low Pacc (~< 10~^) required to 
give the correct comoving densities for A > 2. With this 
in mind, we consider two scenarios which provide the cor- 
rect comoving number density of PBHs (9 x IQ^^Mpcr^ 
assuming a limiting mass of 10^°Mq in section ||), one 
with a reheat temperature of the order = lO^^Gev, 
tilt Us — 1.33 and a second with reheat temperature 
Tr = W^Gev, Us = 1.57 with A = 2 in each case. As 
Fig. H shows most black holes, and in particular those 
that are able to accrete, are then formed immediately af- 
ter reheating, with a slight spread in their masses, which 
we further display in Fig. |lO| . 

We then evolve this initial mass distribution consid- 
ering accretion and evaporation in a quintessence model 

We assume, however, that 



as described in Section [II 



the quintessence field goes off scaling and kinates at a 
temperature of the order T = IMev, since if the field 
continues to scale after this time one ends up with too 
large masses. Nevertheless if the field goes off scaling the 
black holes stop accreting (cf. Eqn ||), and even if the 
field starts scaling subsequently, BHs will not grow sig- 
nificantly again, since their growth pause has rendered 
them subcritical. Hence in what follows what we need is 
simply a brief pause in scaling at a temperature of around 
T = IMev. We shall consider two examples here, one 



' — 

-1.5 

1 — 
0.8 - 




log(M,,(a)/M^„,(a)) 




log(M,„/M„(T,.)) 

FIG. 10. The bottom panel shows the initial black hole 
mass distributions for the two scenarios described in the text, 
with reheat temperatures [Tr) of lO^^Geu (dashed line) and 
lO^Gew (full line). The upper panel shows the same dis- 
tributions comparing them to the critical mass at the time 
of each black hole's formation. A blackhole will accrete if 
{M},h/ Merit) > 1. In each panel the distributions are nor- 
malised to 1 at their peak for ease of comparison. 



with Toff ~ 4Afef, and another with To// — IMev. 

In Figs. 11 and |l2| we find that the projected mass 
spectra at z w 10 (or indeed at any time after Toff but 



before the galactic merging history started). In Fig. |Tl| 
we plot results for PBHs formed at a reheating tempra- 
ture of Tr = 10^ Gev as described above. In Fig. |l^ we 
plot the corresponding distribution for the scenarios in 
which Tr — 10^° Gev. In both cases note that the ex- 
istence of a critical mass for accretion implies that the 
final distribution mimicks the initial one clipped at the 
critical mass. For the first scenario considered the cut off 
is at the peak of the distribution; hence the final distri- 
bution is very skewed. For the second scenario the cut off 
is to the left of the peak - so the final distribution is more 
symmetric. As expected if the field goes of scaling later, 
the final PBH masses are much larger: for Toff — 4 Mev 
we find final PBH masses of the order of 10* — IQ^Mq] 
for Toff — 1 Mev these masses grow to 10^ — 1O^-M0. 

In either case these plots are entirely consistent with 
those obtained from the merging history in section ^ 
(Fig. Quintessence could have therefore provided 

the primordial seeds which then turn into the SMBHs we 
observe today. 
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VI. CONCLUSIONS 
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FIG. 11. The distribution of black hole masses at redshift 
2 ~ 10 for quintessence models deviating from scaling at 
T ~ 5Meu, and IMev (we have used a reheating temperature 
of Tr = 10** Gev and A = 2). In both cases the distribution is 
skewed, reflecting the exitence of a critical mass for growth. 
The later quintessence leaves scaling the larger the BHs' mass. 
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M/Mg 

FIG. 12. The distribution of black hole masses at redshift 
z ~ 10 for quintessence models deviating from scaling at 
T ~ 5Mev, and IMev (we have used a reheating temper- 
ature of Tr = 10^° Gev and A = 2). 



We have demonstrated that a scenario in which pri- 
mordial black holes attain super massive size through 
the accretion of a cosmological scalar field is wholly con- 
sistent with current observational constraints. Such a 
model can generate the correct comoving number den- 
sity and mass distribution for SMBH's, given a standard 
prescription for late-time merging and matter accretion 
and with reasonable choices of tilt. 

Existing schemes explaining SMBHs require very con- 
trived choices of parameters (e.g. [0). The root of all 
the evil lies in the huge mass discrepancy between cosmic 
scale masses (such as those considered in previous stud- 
ies of PBHs) and those of their supermassive cousins - 
the cosmological counterpart of the discrepancy between 
solar mass objects and SMBHs found in astrophysical 
schemes. We attempt to bridge the mass scale gap in our 
model by allowing primordial holes to grow by accreting 
quintessence. Still, we find that we have to switch off this 
process at a carefully tuned time, so that there is not too 
much growth. 

The fact that this requirement is equivalent to re- 
questing that quintessence goes off scaling just before 
nucleosynthesis is of some consolation: such a feature 
is already present in non-minimally coupled quintessence 
models such as those studied by In these scenar- 

ios the present acceleration of the universe results from 
a coupling between quintessence and dark matter. It 
switches on close to the radiation to matter transition, 
but is affected by a long transient. This explains why the 
universe did not start accelerating until nowadays. The 
fact that this transient is symmetric around equality, and 
that equality is an equal number of expansion times from 
us and from nucleosynthesis, then makes the field kinate 
away at nucleosynthesis time. 

In addition, in order to obtain the correct comoving 
density we have to tune carefully the value of the scalar 
tilt, Hs, on the scale of the horizon at the time of pri- 
mordial black hole formation. However this fine tuning 
is a problem with any theory employing primordial black 
holes for astrophysical purposes, such as as candidates 
for dark matter [T6| , and is no better or worse in our 
theory. 

In spite of these fine tuning problems we believe that 
this is an interesting scenario which deserves further 
work. The effects of angular momentum upon the whole 
picture are perhaps the most important next issue to con- 
sider. 
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VII. APPENDIX I - PBH CREATION AND 
QUINTESSENCE 

In this appendix we the examine the effects of 
quintessence on PBH formation using the spherical 
model. The idea is to follow non-linear collapse of a 
super- Jeans size spherical region by modclhng the over- 
dense region as a portion of a Friedmann closed model 
pasted onto a flat model. Setting up this model entails 



studying the dynamics of quintessence in closed mod- 
els. For completeness we shall also consider open models, 
which may be of relevance for modelling the void struc- 
ture of out universe. 



A. Quintessence in open and closed models 



The relevant equations are: 



p + 3-(p + p) 
a 



K 



'-^3-0 = Ae" 
a 



(21) 
(22) 



where dots represent derivatives with respect to proper 
time. It is well know that ii K = Q with p = wp we have: 



(23) 
(24) 
(25) 
(26) 



a oc t^oT^ 
1 

P'^ a3(i+«,) 



n 







2 t 



(27) 



We call this solution the scaling solution. With K ^ 
this is also the solution at early times, when curvature is 
subdominant. 

Without quintessence open models {K = —1) at late 
times become vacuum dominated - the so called Milne 
Universe for which a oc t. However we find that in 
the presence of quintessence the onset of negative cur- 
vature domination leads to another scaling solution: one 
in which curvature and quintessence remain proportional. 
We find that: 



a oc t 

1 

P0 oc — 



-^ = -3 



Pc 

2 



2 

t 



(28) 
(29) 

(30) 

(31) 

(32) 



which implies that open universes at late times are devoid 
of normal matter, but not of quintessence. This type 
of behavior can be understood from the K = scaling 
solution, considering that open curvature behaves like a 
fluid with w — —1/3. In the same way that quintessence 
locks on to matter or radiation in a scaling solution, it 
locks on to open curvature. 
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As is well known, without quintessence closed models 
(K = 1) expand and eventually turn around and collapse 
in a Big Crunch. This type of behavior is not changed 
by the presence of quintessence; however the scaling be- 
havior of quintessence itself is drastically changed. As 
the universe comes to an halt the friction term in the (jj 
equations (that is the term 3^0) is withdrawn and the 
field becomes kinetic energy dominated (i.e. it kinates). 
Hence as the universe turns around quintessence becomes 
subdominant, as it scales like oc in contrast with 
p oc 

However as the universe enters the contracting phase 
what used to be a friction term starts to drive the field, 
since now f < 0. This leading to runaway kination, since 
it is precisely the balance of braking and the slope of 
V{(j)) what usually moderates the balance of kinetic and 
potential energy. As before, kinetic energy domination 
implies = 1 and a faster decay rate with expansion 
p oc However contraction reverses the argument, 

and whatever has stronger dilution rate during expan- 
sion, will have higher compression rate during contrac- 
tion. Therefore, at late stages of collapse quintessence 
dominates. Since curvature and background matter can 
be ignored we have the solution: 



a oc (tc - ty^^ 


(33) 


p.ocl 


(34) 


W0 = 1 


(35) 




(36) 




(37) 



in which is the crunch time. 



VIII. APPENDIX II - QUINTESSENCE AROUND 
A BLACK HOLE 

Given the large Jeans mass of quintessence, its fluc- 
tuations should be very small even in the interior of 
non-linear objects, such as the Solar system. Indeed 
the relativistic restoration force induced by the gradi- 
ents should ensure that quintessence fluctuations remain 
linear, even under highly non- linear gravitational forces. 
One obvious exception is the vicinity of a BHs horizon, 
where </> has to change drastically. 

The interaction of BHs and quintessence then becomes 
a problem of boundary matching: on one side the homo- 
geneous cosmological solution; on the other the infinite 
radius limit of the Schwarzchild solution. This matter 
has been examined in the literature in the context of the 
gravitational memory problem in Brans-Dicke theories 

[38 

We assume that the cosmological scalar field generates 
a much weaker gravitational field than the black hole, so 
that we can impose a Schwarzchild metric, with a quasi- 
stationary mass parameter. Hence the equation for the 
scalar field is the vicinity of the BH is: 



Our analysis will closed mimic that of Jacobson 



1 



A 



dV 



(38) 



where A{r) = 1 — 2M/r, and dots and dashes are deriva- 
tives with respect to time and r respectively. For free 
waves {V = 0) the system separates giving: 



*i?(r) 



(39) 



with: 



B. Implications for structure formation and PBH 
formation 

Qualitatively these results indicate that quintessence 
has a leading role in the strongly non-linear stages 
of structure formation. Voids should be filled with 
quintessence, judging from what happens to the K = —1 
case. Also it would appear that black hole formation 
would be led by quintessence and not by matter (as im- 
plied by the K = 1 case). Interestingly, quintessence ap- 
pears to behave like a stiff fluid (p^ = p^) during collapse, 
so we should be able to simplify the calculations. The 
black hole mass seems to be dominated by the amount of 
quintessence accreted. Also because oc p^a^ oc 
there must be a mass enhancement during collapse (due 
to gravity acting against quintessence's pressure). 

All of these effects may at most introduce factors of or- 
der 1 in the calculations in the main body of the text, and 
so interesting as these results might be we have relegated 
them to Appendix. 



:{r^AR' 



A 



(40) 



or introducing a Kruskal coordinate r*, such that Adr* = 
dr (or r* = r + 2M\og{^ - 1)): 



fR_ 

dr*^ 



2A dR 
r dr 



—+uj^R = 



It is clear that far away the solutions are: 



and that near the horizon they become: 



(41) 



(42) 



(43) 



Focusing on solutions regular on the horizon 7i+ , we in- 
troduce the advanced time coordinate v = t + r*, so that 
the relevant oscillatory solutions take the form e*'^" . 

In addition there are non-oscillatory solutions regular 
at r = 2M, such as those studied by Jacobson. These 
can be obtained noting that (|3^) can be rewritten as: 
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20,™ + + ^ir^A^,r)^r = -V'i<j)) (44) 

For a general y ((/>) separability is lost and numerical work 
is necessary. In order to further the analytical approach 
we model the rolling potential by its gradient at a point 
V = —^(j) (we could add a constant here without loss of 
generahty). Equation ( ^4| ) then becomes: 



2 1,2. ^ 

20,™ + -4>,v + -^[r A4>^r),r = -M 

Setting (f> = f{v) + g{r) leads to / = Bv + D, and 

„2 (J 



9,r 



-fir 



3(r-2M) r-2M r{r-2M) 



(45) 



(46) 



For the solution to be regular at the horizon we must 
have: 



C = AM^ B 



(47) 



Integrating finally leads to: 



B\v 



H r 



r - 2M log r] 

2 



2Afr + 4Af^logr +D 



(48) 



which generalized Jacobson's solution. The general so- 
lution is a superposition of free oscillatory solutions and 
this solution. 

We must now impose a boundary condition of the form: 



= 0c + 



(49) 



where 0c refers to the quintessence cosmological solution. 
If we fix B and D so that: 



2Af log 

^ 2M 



3 V 2 



2A-/r + 4A'/2log — I 



r 

2M J 



we have asymptotically: 



+ 0ct- 



M 2 

— r 
6 



(50) 



(51) 



A small error is made in matching the two conditions, 
but in a quintessence scenario with 



dV 



(52) 



T = T,^ = 4>l (53) 
we finally find an equation for the BH mass: 

47r(2Af)202 (54) 



dM 
IT 



which is the same formula derived by Jacobson. Hence 
we conclude that the potential V{(j)) only has indirect 
impact upon the BH mass, via its effect upon the time 
evolution of the cosmological solution 0c . 



(where A has units of 1/0 and is expressed in Planck units 
and p is the cosmological density) we find thatthe scale 
on which the error becomes significant is the cosmological 
horizon scale. Hence this small error can be neglected. 

Computed the flux of energy through the BH horizon 
associated with this solution: 
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